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Abstract 

The observed gas dynamics in the Milky Way can only be explained by 
I a bar in the galactic center. Such a bar is directly visible in the near- 

. IR maps of the bulge, where it causes a distinctive asymmetric light 

' distribution pattern. Another large-scale structure is the grand-design 

4-arm spiral pattern, most clearly observed in the spatial distribution of 
fH . molecular gas and Hll-regions. Since this spiral arm pattern is regular, 

' it must have existed for at least a few rotations. Traced by molecular 

, clouds, the spiral arms appear strongest in the so-called molecular ring 

• between 4 cmd 7kpc. 

I In order to model the observed gas flow structure, we constructed 

, a model for the stellar mass distribution. For the inner 5 kpc we 

' used the 3D deprojected near-IR light distribution, as observed by the 

^ , COBE/DIRBE experiment, and added an analytical disk model outside 

the box as well as a halo model. With this frozen mass distribution, 
we computed the stationary gas flow for various deprojection parame- 
ters and pattern speeds. For all reasonable parameter choices, we obtain 
' a 4-a.rmed spiral pattern, which can be matched to the observed spiral 

' arms. The model spiral arms, are caused by the forcing of the bar, and 

. two additional mass concentrations at about 4 kpc's on the minor axis 

of the bar, which are likely to be incorrectly deprojected stellar spiral 
arms. 

In the bar region, our model can explain the non-circular motion vis- 
O ' ible in the terminal velocity curve as well as some part of the forbidden 

^ , velocities. Inside the corotation, we also find 4 spiral arms, the near- 

. est arm corresponds to the 3-kpc-arm, although only qualitatively. The 

' m,issing southern 3-kpc-arm at the far end of the galaxy is explained by 

^ , running parallel to another arm. Close to the center, we find gas on cir- 

! ' cular orbits forming a disk. Such a disk has been observed in emission 

K> I of the CS molecule, however only part of the disk appears to be occupied 

^ . by dense enough gas to be traced by CS. 
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1 Introduction 

The Milky Way is our local laboratory of star formation and ISM physics. 
Some observations even have only been possible in the Milky Way, e.g. tur- 
bulent motion and magnetic field in ISM clouds, or search for dark matter 
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candidates. Also, many observable quantities are directly linked to the his- 
tory and structure of the Milky Way, e.g. metalicity gradients (Friedli 1999), 
and distribution of OH/IR stars (see below). A better understanding of the 
detailed mass distribution and gas flow will therefore beneflt many other stud- 
ies. 

Morphology and type of the Milky Way are hard to recognize because 
of dust obscuration and the location of the sun within the disk. Spiral arm 
tangents can be identified in optical and radio brightness distributions (So- 
fue 1973), as well as in molecular gas (e.g.. Dame et al. 1987), HII regions 
(Georgelin & Georgelin 1976), and other tracers; see Vallee (1995) for a com- 
plete list. The 3-kpc-arm, is unusual in this respect: it is very bright in CO, 
but not traced by HII regions. It is therefore possible that it does not form 
stars at the present time. Maybe this is because it has an enhanced turbulent 
velocity or a higher differential shear (Rohlfs & Kreitschmann 1987). The 
3-kpc arm and other peculiar gas dynamics in the galactic center are evidence 
for a bar in the bulge region. The other spiral arms outside the bar region 
within the so-called molecular ring appear to be on almost circular orbits since 
they do not show large non-circular motion in front of the galactic center. A 
comparison of many publications indicates that the Milky Way most likely 
has 4 spiral arms with a pitch angle of 12° (Vallee 1995). 

Another important aspect is, that the true rotation curve of the Milky Way 
can be measured more precisely by a combination of a stellar mass model and 
a hydrodynamical model. When the rotation curve is determined from the 
gas dynamics alone, an axisymmetric model yields an uncommon sharp peak 
at ~ 0.5 kpc (Clemens 1985). An axisymmetric model of the stellar mass 
distribution, however, implies a mostly fiat rotation curve there (Kent 1992). 
A bar naturally solves this issue by explaining the sharp peak with non-circular 
motion caused by orbits elongated along the bar (Binney et al. 1991). It also 
explains the observed nuclear disk of molecular gas inside ~ 200 pc. 

Finally, having a better model for the mass distribution and rotation curve, 
one also obtains better constraints on the amount of dark matter in the solar 
neighborhood. 

Since a short review cannot give a complete coverage of all Milky Way 
models in the literature, we refer the reader to the following papers for similar 
and alternative models: Lin, Yuan & Shu (1969), Mulder & Liem (1986), 
Amaral & Lepine (1997), Wada et al. (1994), Weiner & Sellwood (1999), and 
Fux (1999a). For reviews about spiral structure see: Wlelen (1974), Toomre 
(1977), Binney & Tremaine (1987), and Berlin & Lin (1996). 

2 Using observations to yield a detailed mass 
model 

An early axisymmetric mass model for the inner Galaxy was constructed by 
Kent, Dame & Fazio (1991), by fitting parametric models to the photometric 
near-IR maps at 2.4 /xm (K-Band) which where obtained by the Spacelab In- 
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frared Telescope (IRT) and have a resolution of 1°. By assuming a constant 
mass-to-light ratio for each component Kent (1992) found various possible 
combinations of bulge, disk, and halo components to fit the observed mass 
distribution and kinematics, as well as the gaseous rotation curve outside the 
bulge region. The bulge region was excluded from the fit, because the gaseous 
rotation curve from Clemens (1985) shows signs of non-circular motion due to 
the presence of a bar. Since the dark halo is only observable through its grav- 
itational interaction, Kent used the dark halo component to compensate for 
the mismatch between gaseous rotation curve and the rotation curve implied 
by the disk and the mass distribution in the inner galaxy. He also defined a 
maximum disk model, which minimizes the amount of dark matter required 
by maximizing the mass-to- light ratio for the disk (he found M/L = 1.3). 
This model is particular interesting, because it gives a lower limit for the 
amount of dark matter in the solar neighborhood. One reason why there is 
large uncertainty about the relative contribution of dark halo and disk, is that 
not enough constraints about the vertical distribution of mass in the galaxy 
are known (Dehnen & Binney 1998). 

A major improvement over the model of Kent was made possible by the 
near-IR maps obtained with the DIRBE experiment on board of the COBE 
satellite. By using a foreground dust screen model, Dwek et al. (1995) found, 
that the bulge fight is best fitted by a triaxial fight distribution. This bar 
appears to be elongated with axis ratios 1 : 0.33 : 0.22, the nearer end at posi- 
tive longitudes and inclined by about 20° . Later, an improved dust correction 
and improved parametric model was obtained by Freudenreich (1998). 

Another line of models comes from a non-parametric deprojection method 
introduced by Binney & Gerhard (1996). No parametric model distribution 
has to be prescribed, but an initial model has to be given. The method can 
also be applied to a part of the galaxy, while the model outside that part stays 
fixed as specified by the initial parametric model. Binney, Gerhard, & Spergel 
(1997) applied the method to the DIRBE data and found, that the result is 
robust against changes in the initial model and can reliably recover the 3D 
light distribution in the bulge. They also found additional light concentrations 
on the minor axis of the bar, which they attribute to spiral arm heads in 
that area. While the deprojection method is not able to recover spiral arm 
structure, tests show that spiral arms could produce clumps on the minor 
axis like the ones found in the DIRBE data. The only free parameters in this 
model are the mass-to-light ratio and the orientation of the bar. See Gerhard 
(1996, 1999) for a comparison of these models and further evidence for the 
bar. 

For further modefing, we expanded the mass model from Binney, Gerhard, 
& Spergel (1997) into a series of spherical harmonics to calculate the grav- 
itational potential. The radial density profile in the galactic center follows 
a power law p a r~^'^^ (Becklin & Neugebauer 1968) and this peak in the 
density is not reproduced by the deprojection because of finite resolution and 
smoothing effects. We therefore corrected the zeroth order spherical harmonic 
(monopole) to refiect the power law. 
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The resulting rotation curve falls off beyond 4.5 kpc, because no dark halo 
has been included yet. By construction, the model should reflect the contri- 
bution of all luminous mass in the galaxy without any assumptions about how 
bulge and disk are build together assuming the same M/L. For some models, 
we changed the monopole to allow for a flat rotation curve. 

3 Terminal rotation curve 

The rotation curve of the Milky Way is not directly observable. In external 
galaxies, the rotation curve can be obtained by an observation of Doppler 
shift along a slit across the center. Using parametric models for disk, bulge 
and halo, the rotation curve can then be decomposed. The parametric model 
for disk and bulge must follow the observed surface brightness, while the dark 
matter halo is only constrained by the fit to the rotation curve. Where density 
wave theory is apphcable, i.e. for tightly wound spirals, further constraints 
are available for the dark matter contribution (Fuchs, Mollenhoff & Heidt 
1998). 

In the Milky Way, the rotation curve can only be inferred indirectly from 
the terminal velocity, i.e. the maximum observed radial velocity within the 
galactic plane at a given longitude. Historically, the terminal rotation curve, 
has been used to infer the radial mass distribution by assuming circular ro- 
tation. Inside the solar circle, the terminal velocity is equal to the circular 
rotation curve minus the motion of the local standard of rest (LSR). The 
motion of the LSR can be inferred from the streaming of stars measured by 
Hipparcos (Feast & Whitelook 1997), or the proper motion of Sgr A* (Backer 
& Stramek 1999). Both methods agree reasonable well, and no m = 1 mode 
seems to be present in the center of the galaxy. 

The high peak in the terminal curve, gas on apparently forbidden ve- 
locities, and the large radial velocity of the 3-kpc-arm indicate non-circular 
motion in the inner galaxy. This led to the idea, that our galaxy is actually 
barred and the presence of a bar implies non-circular motion which cannot 
be corrected for without modeling (Gerhard & Vietri 1986). Fortunately, the 
non-axisymmetric light distribution of the bar can be extracted from pho- 
tometric imaging due to perspective effects. These effects are only strong 
enough if the diameter of the stellar system is comparable to its distance. 
This excludes application of similar models for distant objects in the near 
future. A non-parametric method for the deprojection has been recently de- 
veloped by Binney & Gerhard (1996). Its main advantage over parametric 
methods is, that it allows to recover more structure details and the precise 
radial mass distribution. But an apphcation of this method to the DIRBE 
data yielded controversial results (Binney, Gerhard & Spergel 1997; in the 
following: BGS). The recovered radial and vertical structure of bar and disk 
provides a reahstic view of the inner galaxy. But in addition, mass concen- 
trations on the minor axis have been found in the deprojection, which seem 
to indicate that the inner galaxy has signiflcant spiral structure as well. 
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Figure 1: Terminal curve for Kent (1992, long-dashed line), axisymmetric 
BGS model without halo (short dashed line), and gas dynamics in full BGS 
model without halo (soHd Hne and dotted line with higher resolution). 



In Fig. [l] we compare the axisymmetric part of the BGS model (short 
dashed line) with the axisymmetric model by Kent (1992, long dashed Hne) 
and the ^^CO data from Clemens (1985) and HI data from Burton & Liszt 
(1993). Note that Kents model includes a dark halo and assumes LSR motion 
Vq = 234kms"i while the BGS model was plotted for Vb = 220 kms'^ and 
does not include a halo. Increasing the LSR motion increases the gap between 
model and observed terminal curve at the solar circle {I — 90°), but can be 
compensated for with a dark matter component. 

Both axisymmetric models fail to explain the peak in the rotation curve. 
Between 15 and 40° the BGS model fits the data better, further out the 
contribution of the dark halo is missing which is already included in Kents 
model. When we assume a constant rotation curve beyond 5kpc (= 40°), the 
BGS model also matches the data between 40° and 90°. Closed orbits in the 
full BGS potential match the high peak in the rotation curve (Fig. |^). This 
confirms that the deprojection indeed provides a good description of the bar 
potential. Further out, the combination of bar and mass concentration on 
the minor axis complicates the picture. A hydrodynamical model of the gas 
flow in the BGS model is presented below. When the terminal curve for this 
model is compared to the data (Fig. 0; solid line), the result depends on the 
resolution in the gas model. Hydrodynamical forces tend to depopulate the 
orbits responsible for the peak in the rotation curve due to low resolution. A 
model with higher resolution matches the peak better (Fig. |^; dotted line). 
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Figure 2: The high peak in the observed terminal velocity curve (points) is 
reproduced by the a;i-orbital family in the BGS model. 

Whether bars and spirals are independent dynamical entities is not yet 
clear and may vary from galaxy to galaxy. Numerical A/'-body simulations 
show that a bar can coexist with a spiral of much lower pattern speed (Sell- 
wood & Sparke 1988). On the other hand, detailed modeUing of the observed 
gas dynamics in NGC 1365 was possible using bar and spiral perturbation 
with a single pattern speed (Lindblad, Lindblad & Athanassoula 1996). 

Moreover, the spiral mode is most likely growing and decaying within a 
few rotations of the galaxy (Lin & Shu 1964). For our gas model, we assumed 
that these additional light concentrations can be used as a first approximation 
to the contribution of spiral arms in this region and that the spiral pattern is 
rotating with the same speed as the bar and is stationary. 

A more advanced deprojection technique combined with spiral arm mod- 
eling will be used to shed more light on this issue shortly (Bissantz & Gerhard 
2000). This will also allow us to study models with independently rotating 
bar and spiral mode. 

4 Steady state / Stationary gas flow solution 

If the gas is allowed to relax in the potential of the galaxy, i.e. the mass 
distribution in the galaxy does not change over a few rotational periods, we 
expect the gas to form a stationary flow pattern. Since the gas looses energy 
in collisions, it arranges itself to follow nested closed orbits which have no 
intersections. In barred galaxies, orbits are organized in famihes. The so- 



6 



called Xi-orbits are elongated along the bar and exist from somewhat inside 
the corotation annulus to the center. Inside the bar region there is a special 
cusped xi-orbit, which is the last xi orbit the gas can settle on. Inside this 
orbit, the xi orbits in our potential develop loops with self-intersections. Gas 
on such orbits, would suffer coUisions causing it to leave the orbit. In Hnear 
theory, the gas has to pass through the inner Lindblad resonance (ILR), which 
causes a phase shift of 7r/2 in the orbits. The gas follows this phase shift 
with a delayed response, and switches to the X2 family of orbits, which are 
elongated perpendicular to the bar and exist only inside the ILR. In general, 
further ILRs may exist, but there is no evidence for further ILRs in the BGS 
potential. The transition at the ILR is not smoothly, but is accompanied by 
a shock in the gas flow. In strong bars, the shocks are straight Hues, which 
in real galaxies have been identified with dust lanes and velocity jumps. Due 
to the non-axisymmetric distribution of gas caused by the shocks or spiral 
arms, the bar imposes a gravitational torque on the gas, which causes the gas 
to flow in (out) when the traihng spiral arm is inside (outside) corotation. 
For leading spiral arms, the torque would be reversed. For this reason, gas is 
slowly depleted in the corotation region and accumulates within the ILR to 
form a ring or disk of gas on X2 orbits. 

For hydrodynamics we chose the smooth particle hydrodynamics (SPH) 
method (Benz 1990; Steinmetz & Miiller 1993). In SPH a continuous gas 
distribution is approximated by a spatially smeared out particle distribution. 
The advantage of this method is, that self-gravity can easily be included. 
Further numerical details of this simulation are given in Englmaier & Ger- 
hard (1999). In addition to the free parameters from the deprojection, our 
gas model needs only three additional parameters: the corotation radius or 
equivalently the pattern speed, the effective sound speed, and the LSR motion 
for calculation of (Iv) diagrams. Our best model is shown in Fig. ||. 

The corotation radius in our model is tightly constrained by the observed 
(Iv)-diagram as follows. Since the 3-kpc-arm is clearly in non-circular rotation, 
it must be inside the corotation radius. By comparison, the 4-armed spiral 
pattern in the molecular ring between 4 and 7kpc shows no large deviation 
from circular motion. This observation can only be reproduced in our model, 
if the corotation is between 3 and 4kpc. The angular extend of the 3-kpc-arm 
is best reproduced for corotation at 3.4 kpc. 

The M/L ratio and the LSR motion are fixed by a fit of the terminal 
velocity curve to the model. Actually the circular rotation velocity for the 
LSR is kept constant to 200kms~^ times a scaling constant which depends 
on M/L. The final value for the LSR motion is 208 kms^^ after the best 
model fit has been determined. While this value is lower than the current 
best estimate for the LSR motion, we note that the M/L value and other 
model details do not depend strongly on this parameter. The difference in 
the terminal curve can be compensated for by a dark halo. 

The gas fiow model confirms the allowed range for the bar incHnation 
found by BGS (25 ± 10). Models with 20° or 25° yield a better overall fits 
than models with 15° or 30°. A better model for the stellar spiral arms may 
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Figure 3: Comparison of our model with observed spiral arm tangents (lines, 
averaged values) and HII regions (stars) and giant molecular clouds from 
Dame et al. (1986; circles). Rotation is clock wise, and the model is viewed 
from the north galactic pole. The bar inclination is 20° and corotation at 
3.4 kpc. This model also includes a dark halo. 

improve this situation. 

5 Identification of spiral arms in the (Iv) -diagram 

A full sky survey of emission from atomic hydrogen and molecular species 
Hke ^^CO allow mapping of large scale structure such as spiral arms (Fig. ^). 
The longitude-velocity diagram (Iv-diagram) shows traces of spiral arms as 
crowded regions. By assuming circular rotation, one could map the observed 
(Iv)-position back to real space, however, this leads to serious errors close 
to tangential points (Burton 1971). Nevertheless, with some success, the 
principle spiral arm structure has been inferred from Iv-diagrams. First, this 
was done by Oort, Kerr & Westerhout (1958), later the 4-armed structure 
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Figure 4: The (Iv)-diagram for ^^CO from unpublished data by Dame et 
al. (1998). This figure contains all emission integrated over latitudes between 
b ~ —2° and b — 2°. The grey scale is adjusted such as to emphasize spiral 
arm structures. 

was first claimed by Georgelin & Georgelin (1976) using the (Iv)-diagram for 
large Hll-regions. This picture has later been refined by Caswell & Haynes 
(1987), Lockman (1989), and others. In Fig. || we compare various known 
Hll-regions and giant molecular clouds from the literature with our model. 

HII- regions are a particular good tracer for spiral arms, as is known from 
observations of external galaxies. The 4 spiral arms are also prominent in 
the so-called molecular ring between 4 and 7kpc. This area corresponds to 
a bright band of ^^CO emission between about 30.. 60° on the northern side 
of the galaxy (left in Fig. ^ to —(30. .60)° on the southern side. As Solomon 
et al. (1987) pointed out, the spiral arm tangent at around +30° is actually 
split into two components at 25° and 30°. A summary of available spiral arm 
tangent data is listed in Table 1. 

The molecular emission also shows an arm apparently not traced by HII- 
regions: the 3-kpc-arm. This arm shows large non-circular motion, as it 
passes in front of the galactic center (where it appears in absorption) with 
-1-54 km s~^ radial velocity. All other arms show much less radial velocity in 
front of the galactic center. 

Our gas model qualitatively reproduces many of the prominent features 
in the observed (Iv)-diagram (Fig. The 4 spiral arms are found to be 
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29 
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-32 
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24, 30.5 


49.5 


-50 


-30 
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25, 32 


51 








(4) 


25, 30 


49 








(5) 


24, 30 


47 


-55 


-28 




(6) 


32 


46 


-50 


-35 
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32 


48 


-50,-58 


-32 


-21 


(8) 


29 






-28 


-21 








47 


-Ol 


9n 

-zu 


average 


oU 




-Ol 


-ol 




Model 












(a) 


- 25 


54 


-44 


-33 


-20 


(b) 


- 30 


50 


-46 


-33 


-20 


(c) 


- 29 


51 


-47 


-34 


-22 



Table 1: Measured spiral arm tangents in inner galaxy: (1) atomic hydrogen 
(Weaver 1970; Burton & Shane 1970; Henderson 1977); (2) integrated ^^CO 
(Cohen et al. 1980; Grabelsky et al. 1987); (3) ^^CO clouds (Dame et al. 
1986); (4) warm CO clouds (Solomon et al. 1985); (5) HII- Regions using 
H109-a (Lockman 1979; Downes et al. 1980); (6) ^^Al from massive stars 
(Chen et al. 1996); (7) Radio 408 MHz (Beuermann et al. 1985); (8) 2.4 /xm 
(Hayakawa et al. 1981); (9) 60 ^m (Bloemen et al. 1990). Compared to 
some models: (a) Rc = 3.4 kpc, bar inclination Lp = 20°, without halo; (b) 
Rc = 3.4kpc, (fi = 20°, with halo vq = 200 kms'^; (c) Rc = 3.4kpc, = 25°, 
with halo vq = 200 km s^^. 



embedded in the molecular ring with about the right tangential directions (see 
Table 1). Our model, however, assumes perfect point symmetry in the gas, 
while the galaxy is not. Hence a perfect match cannot be expected in such an 
ideaHzed model. Inside the molecular ring, gas is depleted in the corotation 
region and the spiral arms show a gap (see Fig. Within corotation the 
spiral pattern continues followed by the usual gas flow conflguration in a 
barred galaxy (see previous section). One of the spiral arms is similar to 
the 3-kpc-arm in the galaxy, alas with too small an expansion velocity. In 
this region, the gas flow may be disturbed by the stellar spiral arms. In a 
modifled model where we approximate the spiral arm gravity by the gaseous 
model spiral arms, but spatially smeared out to account for wider stellar spiral 
arms, the 3-kpc-arm can be reproduced quantitatively (Englmaier & Gerhard 
1999). Our model also predicts the location of the arm corresponding to the 3- 
kpc-arm on the far side of the galaxy. It happens to fall close to the molecular 
ring and indeed this arm tangent was found to be split into two components 
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by Solomon et al. (1987) and is also visible at I = 25°, 30° in Fig. A similar 
position for the far 3-kpc-arm was suggested by Sevenster (1999), but see Fux 
(1999a) for an alternative explanation. A nuclear disk with ~ 200 pc radius is 
formed in the model by gas on X2-orbits. Size and rotational velocity of this 
disk depends on the enclosed mass which has been adjusted by the central 
density cusp profile. While the disk matches the observed CS emission (Stark 
et al. 1991), only part of the disk is occupied by dense enough gas to show 
CS emission. 

6 OH/IR stars as a fossil dynamical record 

Sevenster (1997, 1999) made a complete sample of the OH/IR stars in the 
bulge region between -45° < I < 10° and |6| < 3°. The OH/IR stars are 
giant stars that lose matter in the so-called asymptotic giant branch (AGB) 
superwind phase. While these stars are quite old, i.e. ~ 1 to 8 Gyr, they 
quickly go through an evolutionary phase characterized by OH maser emis- 
sion in the IR. The maser emission lasts only for 10^~^ yr which is short 
compared to dynamical timescales. The sample of OH/IR therefore provides a 
snapshot of the dynamically evolved star formation regions ~ 1 to 8 Gyr ago. 
Sevenster (1999) found, that the Iv-diagram for the sample shows a striking 
correlation between the OH/IR stars and the 3-kpc-arm observed in ^^CO. A 
young group of about 100 to 350 Myr old OH/IR stars follows this arm very 
nicely, while the older OH/IR population is more equally distributed. This 
appears to be in contradiction with dynamics, since a circular orbit at 3 kpc 
takes ~ 2nr/v ~ 80 Myr to complete. Any trace of the spatial distribution 
of star formation would have been wiped out. However, if these stars were 
formed in one of the spiral arms which itself rotates, then only the relative 
orbit between stars and arm would matter. For a constant pattern speed 
equal to the bar pattern speed, say Q ^ 60 kms^^ kpc^^, the orbit in the 
rotating frame of the arm would take about ~ 27rr/(w — fir) ^ 0.5 Gyr to 
complete. This much longer timescale may allow stars formed at the same 
place to generate a distribution similar to the observed one. Furthermore, the 
3-kpc-arm is not a strong shock, making it possible for stars formed within to 
drift away rather slowly. 

In order to estimate the dynamical constraints of our model on the ob- 
served OH/IR star distribution, we picked our best model and selected gas 
particles which are in compressed areas, i.e. in spiral arm shocks (see Fig. ^). 
We evolved these particles as test particles in the background potential of the 
model, and plotted snapshots for the particle distribution at later times. It 
turns out that, although OH/IR stars may form a transient arm like feature, 
it will only last for a few 100 Myr. In Fig. ||, we show the Iv-diagram for 
100 Myr old stars. There is a 3-kpc-arm-like distribution of stars, but these 
stars come from the corotation region at the minor axis, not the 3-kpc-arm. 
Just inside corotation, a star formed from the gas moves faster than the bar 
pattern, and this gives raise to a coriolis force pulling the star inwards. A 
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Figure 5: Location of star forming regions of compressed gas in face-on view 
(upper) and seen in the {lv)-diagram (lower). For each particle a 2Myr long 
trace of its orbit is shown, which is approximately the duration of the OH/IR 
phase. The bar is aligned with the diagonal Hne from the lower left to the 
upper right (compare with figure ||) . The sun to galactic center Hne is incHned 
25° to the bar. The location of the sun is just outside the diagram at a; « 
— 2.7kpc, y « —7.5 kpc. The 3-kpc-arm is the horizontal trace of stars starting 
at the lower left end of the bar. For the (Iv)-diagram a circular velocity 
of 220kms~^ has been used for the transformation in the LSR. The 3-kpc- 
arm passes the galactic center {I — 0°) with the velocity ~ — 40kms~^ and 
continues to I — —21°, v = — lOOkms^^. 
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Figure 6: Same as figure ||, but 100 Myr later. Stars created in the 3-kpc-arm 
have formed elliptical rings ahgned with the bar which do not follow the 3- 
kpc-arm in the (Iv)-diagram, but are restricted to |Z| < 10°. However, a trace 
of stars originally created in the corotation region forms a bridge between 
corotation and ends of the bar maintaining a 3-kpc-arm alike structure in the 
(Iv)-diagram. 
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group of freshly formed but isolated stars would form a stream which may be 
coincident with the 3-kpc-arm. As we discussed above, the 3-kpc-arm does 
not contain Hll-regions and presumably, it docs not form stars because it is 
dynamically hot. Stars which are formed in the 3-kpc-arm are quickly dis- 
tributed on closed orbits, it is therefore possible that OH/IR stars are indeed 
distributed in rings indicating radii of more compressed gas in our model, and 
these orbits, which are not circular orbits, may be indirectly correlated with 
the 3-kpc-arm, because arms which are driven by a bar, often start with a 
local density enhancement and strong star formation on the bar major axis. 
However, we would expect those rings to be inchned with respect to the 3- 
kpc-arm, unless the 3-kpc-arm is not stationary. 

A similar result was found by Fux (1999b), where the observed OH/IR star 
distribution could be reproduced for about 40 Myr old stars. The OH/IR stars 
in his model originate in the far end of the bar in an area of enhanced density. 
Contrary to our model, Fux's model has the additional advantage, that the 
3-kpc-arm changes rapidly and thus the gas is on more ballistic orbits. He 
did not show any evidence, however, that the coincidence of stars and gas can 
be maintained over multiple rotations as required by the observations. While 
both models are optimized to fit various aspects of the Milky Way, they both 
suggest that OH/IR stars are not directly tracing spiral arms. 

We therefore consider it more likely that the stars formed close to coro- 
tation where timescales are comparable to the age of these stars. It will be 
interesting if further evidence can be found, to show that OH/IR stars follow 
closed orbits and whether these favored orbits have changed with time. 

7 Discussion 

We presented a coherent description for the structure and gas dynamics in 
the Milky Way reaching from small scale ~ 100 pc to large scale ^ kpc. Our 
model provides a hnk between photometric and gas dynamical observations, 
but does not intend to be a a self-consistent treatment of formation and 
evolution of the Milky Way. Idealized assumptions about symmetry and in 
the description of the ISM have been made to isolate generic from peculiar 
features. 

We find, that the 3-kpc-arm can be explained qualitatively by the forcing 
of the bar. The bar may also be responsible for some of the spiral structure 
observed in the gas. Strong bars can drive spiral arms to and somewhat 
beyond the outer Lindblad resonance (Mulder 1986). The deprojection of the 
near-IR bulge light distribution introduced additional structure on the minor 
axis of the bar at around corotation. When this structure is included into the 
model, four spiral arms are driven in the gas similar to the observed spiral 
arm structure in this region. The pattern is very similar to a model by Mulder 
& Liem (1986), but the location of the sun is different. Interestingly, Mulder 
(1986) mentions that this pattern can be driven by a strong bar or by a weak 
bar plus spiral arms beyond corotation. 
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When the gravity of the spiral arms is included in the model, the 3-kpc- 
arm can also be explained quantitatively. This finding, and in addition our 
inability to find a fitting model for the 3-kpc-arm in models without the minor 
axis structure demonstrates, that the inner galaxy is affected by both bar and 
spiral arms. 

While spiral arms in barred galaxies usually start at the end of the bar, 
they do not need to be driven by the bar (Sellwood & Sparke 1988), and might 
be dynamically independent with much lower pattern speed. No steady state 
but maybe a periodic steady state can be found in the gas flow of such a 
model. Independent stellar spiral arms further comphcate the gas dynamics 
to some extent, because they cause non-linear velocity jumps in the gas and 
will not just add to the structure formed by the bar. From density wave 
theory (Lin & Shu 1964) we know, that spiral arms are not stationary but 
form and decay in a few rotational periods. With a lower pattern speed, the 
spiral pattern would extend to larger radii. A model for the galaxy, where 
spiral arms and bar have different pattern speeds has still to be constructed, 
yet. 

Considering all these objections, it is even more surprising how well our 
model explains the spiral structure. Both the bar and the spiral pattern 
impose similar constraints on the orientation of the pattern, i.e. the phase 
angle. Nevertheless, this overlap might be a mere coincidence and a better 
treatment of bar and spiral mode is underway to investigate other possibilities 
(Bissantz et al, in prep.). 
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